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Abstract 

We investigate near the point of glass transition the expansion of the free energy corresponding 
to the generalized Sherrington-Kirkpatrick model with arbitrary diagonal operators U standing 
instead of Ising spins. We focus on the case when U is an operator with broken reflection symmetry. 
Such a consideration is important for understanding the behavior of spin-glass-like phases in a 
number of real physical systems, mainly in orientational glasses in mixed molecular crystals which 
present just the case. We build explicitly a full replica symmetry breaking (FRSB) solution of the 
equations for the orientational glass order parameters when the non-symmetric part of U is small. 
This particular result presents a counterexample in the context of usually adopted conjecture of 
the absence of FRSB solution in systems with no reflection symmetry. 
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I. INTRODUCTION 



The theory of spin glasses has been developed as an attempt to describe unordered 
equilibrium freezing of spins in actual dilute magnetic systems with disorder and frustration. 
This problem was soon solved at the mean- field level [see also Ref. |J for a review]. 

Below the Almeida - Thouless line [[| the replica symmetric solution was shown to be 
incorrect. Parisi proposed the method of replica symmetry breaking (RSB) step by step 
with the limit — full RSB (FRSB) when glass order parameter becomes a continuous non- 
decreasing function q(x) of a parameter < x < 1. It provides the hierarchical distribution 
of pure states overlaps probability P(q) through P(q) = dx/dq. This approach allows to 
describe the main features of the experiments on spin glasses. 

A number of generalizations of the Sherrington and Kirkpatrick (SK) model with Ising 
spins have been considered. But still the problem how RSB and FRSB solutions do occur 
remains relevant and far from being completely understood f7HT3|. It was shown that the 
violation of the replica symmetry is correlated with the symmetry properties of the Hamil- 
tonian |I4|^i9|. There is a conjecture that in the absence of the reflection symmetry it is not 



possible to construct a continuous non-decreasing function q(x) and so, the FRSB solution 
does not exist in this case, or at least does not occur in the point of RS solution instability. It 
is possible that in some of these models different stages of replica symmetry breaking occur 



not in the first symmetry breaking point as in the SK model (see, for example, Refs. [£(J, [21 
where the Potts model is considered). 

In the literature the absence of reflection (or time reversal) symmetry usually was incor- 
porated in the structure of the Hamiltonian. It is so, for example, for Potts spin glasses and 
for p-spin spin glass with the interaction of p Ising spins. However, there is another way 
to break reflection symmetry - that is to incorporate the breaking in the character of the 
operators themselves while the structure of the Hamiltonian remains two-particle as in the 
SK model. As far as we know the present paper is the first one to deal with RSB in such a 
situation. 

It is worth to notice that, while, say, p-spin model has a physical meaning only as a 
prototype of structural glasses, there exists a number of real physical systems which can be 



described in terms of two-operator random systems [17, 22, ^3] of the SK type where the 



absence of reflection symmetry is caused by the characteristics of the operators U themselves. 



Let us list some physical examples of such systems. They are just described by the 
Hamiltonians of the SK type but with Ising spins changed for non reflection symmetrical 
operators U. The investigation of replica symmetry breaking presented in this paper can 
give some new information about low-temperature behavior of these systems. 



For example ortho — para-hydrogen mixed crystals and in Ar — N 2 present mixtures 
of spherically symmetric molecules and momentum-bearing molecules. The corresponding 
orientational quadrupolar glass was investigated on the base of the S-K type Hamiltonian 
with U = Q, where Q = 3J| — 2, J = 1 [23 1. Another example of a SG-like phase 
in molecular crystal is presented by pure para — H2 (or ortho — D?) under pressure. The 
possibility of orientational order in systems of initially spherically symmetric molecule states 
is due to the involving of higher order orbital moments J = 2,4... in the physics under 
pressure. The frustration and disorder give the basis to the investigation of quadrupole 
glass with J = 2. Now U = (1/3) (3 J 2 Z - 6). Such a theory was constructed in Ref. |22| . 
Another interesting example of spin-glass-like phase is presented by the orientational glass 
that appears in molecular Cqq at 90K. Now the sample is not a mixture and the form of the 
molecules is almost spherical. Nevertheless the existing anisotropy of the potential causes a 
frustration. The dependence of the anisotropic part of the potential on mutual orientation 
of molecules has two pronounced minima and so the effect of a mixture is obtained. This is 
the base of the theory with U being a superposition of cubical harmonics p5f . 



The aim of this work - to find out for which class of ^/-operators that do not have the 
reflection symmetry the problem can be described by the Parisi FRSB scheme. In this 
paper, we do the first step in this direction. In Ref. |26| it was shown that in the case of 



arbitrary ^/-operators satisfying the condition of the reflection symmetry there is full replica 
symmetry breaking. Keeping in mind that the FRSB solution is valid for operators with 
the reflection symmetry we seek for the FRSB solution for operators U represented by the 
sum of a symmetric operator and a small perturbation that has no reflection symmetry. 
This investigation is close in spirit to the problem of building the FRSB solution for the SK 
model in a weak field HI. 



II. THE MODEL 



We start with the Hamiltonian 



H=-\Y. J V (j rU J 



(1) 



where arbitrary diagonal operators U are located on the lattice sites i. The quenched 
interactions J,j are distributed with the Gaussian probability 

{J hJ fN 
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where N is the number of sites. Using replica approach we can write in standard way the 
free energy averaged over disorder in the form 
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Here 
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where t = J/kT. The saddle point conditions were used to define the glass order parameter 



q aP = Tr 



and the auxiliary order parameter 
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In the RS approximation from the extremum condition for the free energy for the glass 
order parameter q RS we have: 
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Here 



and 



0rs = ztyfq^U + f 
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72. ' V 2 / (10) 
The stability of the RS solution can be tested by the investigation of the gaussian fluc- 
tuations contribution to the free energy near this solution. The solution is stable while all 
the eigen modes of the fluctuation propagator are positive. The most important mode is the 
so-called replicone mode 0, 2?| since only its sign is usually sensitive to the replica symme- 
try breaking degree and to the temperature. For example, the replica symmetric solution 
is stable unless the corresponding replicon mode energy A( RS ) rcp i > 0. The RS-solution can 
break at the temperature T c determined by the equation A( RS ) rcp i = 0, where 
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The equation A( RS ) rep i = is nothing else that the bifurcation condition for 1RSB equation 



for glass order parameter |23], 2J|, i.e., as the condition that a small solution with 1RSB can 
appear. Analogously the other \( n -i)RSB and the bifurcation condition at the n-th stages of 
the replica symmetry breaking are related. 

We break the RS once again and obtain the corresponding expressions for the free en- 
ergy and the order parameters. The bifurcation condition A( 1RSB ) rC pi = determining the 
temperature follows from the condition that a nontrivial small solution for the 2RSB glass 



order parameter appears. [23|, |27, 29 1 We have: 
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where 6 1RSB is the analog of (|]) for the first stage of RSB and m - 1RSB order parameter. 



III. FRSB SOLUTION 



Let us consider first a generalized model defined by the Hamiltonian ([!]) with reflection 
symmetrical operators U. The reflection symmetry implies that for any integer k, 



Tr 



0. 



(13) 
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In the RS-approximation we find the solution g RS that is zero at high temperature. The 
bifurcation condition in this case is: 



l-#«L(*c) = 0. (14) 

This equation coincides with A( RS ) rcp i = [see, e.g., Ref. ||. It is zero high temperature 
solution that bifurcates. At T < T c certain nontrivial 1RSB solutions appear but they are 
unstable. 

Investigating 1RSB, 2RSB, 3RSB, nRSB, and so on, we see that the equations for the 
glass order parameters always contain the quantity 

Tr[U exp(0 nRSB )]/Tr[exp(6> nRSB )]. (15) 

Here 9 nnSB are the analogs of M) for higher stages of RSB (see Ref. |SB]for details). Therefore, 
one of the solutions of this equation is trivial at each RSB-step, and the appearance of the 
nRSB solution can be regarded as the bifurcation of the trivial (n — 1)RSB solution. In 
this case, the equation A nRSB = coincides with the corresponding branching condition 
(fUj). This means that in any case, the nRSB solutions at different stages of the symmetry 
breaking can exist at temperature T < T c determined by this bifurcation condition, and so 
we always can look for FRSB solution. Writing the free energy as a series over 5q al3 near 
T c (up to the fourth order) we obtain q{x) = cx in the leading approximation [a similar 
procedure was described in details in Ref. [2~6|j . 

If the operators U do not have the reflection symmetry, Tr ?/( 2fc+1 ) ^ 0, then the glass 
freezing scenario is different from the previous case. The characteristic properties of system 
develop themselves already in the replica symmetry approximation. The nonlinear integral 
equation for the RS-glass order parameter simply has no trivial solutions at any temper- 
ature because the integrand is nonsymmetric due to the cubic terms in the free-energy 



expansion P7|| . There is a smooth increase in the RS order parameters as the tempera- 
ture decreases. The bifurcation condition A( RS ) rep i = (|Xl|) defines the point T c where the 
RS-solution becomes unstable. 

Similarly to the previous case, considering respectively 1RSB, 2RSB, 3RSB, nRSB, 
we find that the equation A nRSB = always has the solution which determines the point T c 



and coincides with the solution of equation A( RS ) rcp i = [^3], ^ EH- ^ ^ s important that it 
is the non-zero solution that bifurcates. 
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Looking now in general at the free-energy series over the glass order parameter we see 
that the series contain explicitly the terms which can be classified by the reflection symmetry. 

To estimate the form of the FRSB-solution near the bifurcation point T c at which it ceases 
to coincide with the RS-solution (i.e., in neighborhood of T c ), we expand the expression for 
the free energy up to the fourth order (inclusively), assuming that the deviations 5q a ' 3 from 
q KS are small. We believe that one can neglect the changes of the order parameter w KS (see 



Ref. |29] where it is directly shown for 1RSB). We obtain the deviation AF of the free energy 
from its RS part: 
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a,/3,7,<5,/x a,/3,7,5,/i J 

where i = t c + At. The prime on the sum means that only the superscripts belonging to the 
same <5g are necessarily different in ^ . The expressions for coefficients W, L, D are given 
in Appendix. All coefficients depend only on RS-solution at T c . Note that l—t 2 W = A( RS ) r epi- 
The expression (|TjJ) includes a part without the reflection symmetry, namely, the terms with 
odd number of identical replica indices (see also JIT|--|i"5[ ) . 

To write the free energy as a functional of q(x) we use the standard formalized algebra 
rules [§, [|. The properties of this algebra were formulated by Parisi for Ising spin glasses. 
In our case, the expansion of the generalized expression for the free energy Eq. (|i"6|) includes 
some terms of non-standard form. Those terms are not formally described by the Parisi rules, 
but can be easily reduced to the standard form. To do this, we compared the corresponding 
expression, consistently producing 1RSB, 2RSB, ... symmetry breaking. For example, 

lim - Y ' 5q aP 5q ai $q aS = lim - V ' 5q afi 5q^5q ps (17) 

n-*0 71 ^ — ' n->0 fl *■ — ' 

a,/3,7,<5 q,/3,7,5 
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The equation for order parameter follows from the stationarity condition (8/8q(x))AF = 
applied to the free energy functional. The resulting complicated integral stationarity 
equation can be simplified using the differential operator O = ^J^-r Jj:, where q' = ^j^-- As 



d.v 



a result, we obtain: 



t 6 {B 4 -B 3 x} + t 8 l-D 46 x(q)+ 
D 47 



—4q(x)x — Ax(q) + Ax / dyq(y) 

o 



D 



33 



Aq(x)x + 2(g) - 2 / dyq(y) 
o 



AD 2 q{x) } = 0, (18) 



where (q) = f Q dyq(y). Our results agree with those obtained in Ref. ^ in the case when U 
are Ising spins and {U) 2 = 1. 

Differentiating the equation Eq. (|l^) we obtain in the leading approximation: 



tH [-D 2 + D 33 x - D 47 x 2 



(19) 



In deriving the equation Eq.([T9]) we have neglected in the numerator of the members 
terms of the form q(x) ~ r = (T c — T)/T c compared with the constant B 3 . This is true 
for the consideration of Parisi {§-||. Since the model with the operators U with reflective 



symmetry ( i.e. Tr 



0) exactly the same as the model of the Ising spins |26] , for 



reasons of continuity, we consider operators who have a very small Tr 
It is easy to see that 
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It follows from the Cauchy-Schwarz inequality that the expression ^Tr?7 2 e^ RS j ^Tre^ RS j > 
TiUe^Y follows from (£ n A„ 2 ) (E„^ n 2 ) > (£ n A nJ B n ) 2 - 



We have also 
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If U are operators with reflection symmetry, then D 33 = 0, and the denominator of 
is positive. So, in the case of reflection symmetry q' > 0. 

If the operators U have no reflection symmetry, then .D33 7^ and the denominator of ( |18|) 
can be negative. However, for reasons of continuity, one can imagine that it is not always the 
case. As an example, let us consider the operators U = S + t]Q where 77 is small. Here S is 
the ^-component of spin (for S = 1) taking values (0, 1, —1). While Q is the axial quadrupole 
moment, Q = 3S 2 — 2, and it takes values (—2, 1, 1) (see, e.g. [^6|). The operators Q and S 
have the following properties: Q 2 = 2 — Q, 3S 2 = 2 + Q, and QS = SQ = S. So the algebra 
of these operators is closed. The operator S has the reflection symmetry while Q has not. 
FRSB is valid for arbitrary reflection symmetric operators p6[ , in particular, for 5". Let 
us note that the operator v^S* = V is a second component of the quadrupole momentum 
operator V = S 2 — S 2 considered in the problem of anisotropic quadrupolar glass. 

The detailed calculation was performed and it was shown that for small 77 there is the 
vicinity of Tc where q'(x) remains to be positive. The area where q(x) depends on x is small 
and the function q(x) is small ~ r = (T c — T)/T c . 

Thus we can construct correct FRSB solution with q' > in the absence of reflection 
symmetry of operators U. 

It is a possibility that at lower temperature the -D33 term occurs to be larger so that the 
derivative ([18]) becomes negative and the FRSB ceases to exist. This fact can give rise to a 



"reverse" behavior of the system as compared with "standard" case proposed in Ref. fJU] for 
Potts glass models. 
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IV. CONCLUSIONS 



So, we have considered a model with two-particle interaction where the absence of reflec- 
tion symmetry is caused by the characteristics of the operators U themselves. FRSB is first 
described in such a system. An expansion for the free energy of our generalized SK model 
with arbitrary operators U standing instead of Ising spins is investigated near the RSB tran- 
sition point T c . The principal prescription for obtaining a full replica symmetry breaking 
solution is derived in general case. In a case when U is a reflection symmetric operator with 
a nonsymmetric perturbative part the FRSB solution is constructed explicitly. 
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VI. APPENDIX 



W =(U{Uj) - 2(UfU 2 U 3 ) + (UiUzUzUt) = / dz 



Tr [ U 2 e § ™ 



Tr e *s 



Tr Ue § ™ 
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Tr e^ R s 
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Notation used below are obvious from the equation (|23|). The coefficients of the terms of 
the third and the fourth order we write out only those that are included in the final equation. 



L={U?U 2 U 3 )-(UiU 2 U 3 U 4 ). 



(24) 



B 4 = ^{UiWsUaUA) - (UfU 2 U 3 Ujl 6 ) + ^ 1 3 ^ 2 f> 3 f>4) + 

\{UlUtW±) - \{UlU 2 A) + ^(UlUlYi (25) 
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D 31 = -J(tt) + ^(^#4> + ^ 

^{ulutuluA) - + 7(u?u$u 3 u A u 5 u (i ) + l(u?u 2 u 3 U4U 5 u 6 )- 

7(UlU 2 U 3 U 4 U 5 U 6 U 7 ) + 2(U 1 U 2 U 3 U A U 5 U 6 U 7 U 8 ) ] (26) 

D 33 = -l{UlUlUl)+l{^ 

21 

(U?U 2 U 3 U 4 U 5 U G ) + —(UlUlU 3 U,U 5 U 6 ) - A(U!U 2 U 3 U 4 U 5 U e U 7 ) + (UMUsUiUsUeUrUs); 

(27) 

l(U!U 2 U 3 UiU 5 U 6 ) + liUlWzUAUsUr) ~ (uAu^U.UeUrUs). (28) 
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